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We have introduced a class ol exactly soluble Hamiltonian with either SO(2n + 1) or SU(2) 
symmetry, whose ground states are the SO(2n + 1) symmetric matrix product states. The hidden 
topological order in these states can be fully identified and characterized by a set of nonlocal string 
order parameters. The Hamiltonian possesses a hidden (Z2 X Z2) 71 topological symmetry. The 
breaking of this hidden symmetry leads to 4™ degenerate ground states with disentangled edge 
states in an open chain system. Such matrix product states can be regarded as cluster states, 
applicable to measurement-based quantum computation. 
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Quantum spin systems have shown many fascinating 
phenomena and stimulated great interest in the past 
decades. Based on semiclassical argument, Haldane pre- 
dicted that there is a finite excitation gap in the ground 
state of an integer antiferromagnetic Heisenberg spin 
chain |l| . This intriguing feature of quantum spin chains 
results from the breaking of a hidden topological sym- 
metry embedded in the valence bond solid state pro- 
posed by Affleck, Kennedy, Lieb, and Tasaki (AKLT) 
2]. The valence bond solid is a matrix product state 
in one dimension. It shows a striking analogy to the 
Laughlin ground state for the fractional quantum Hall 
effect [3j, |j| . To characterize this topological symmetry, a 
set of nonlocal string order parameters were introduced 
[a, ly] . These string order parameters provide a faithful 
quantification of the hidden antiferromagnetic order of 
the S = 1 Heisenberg model. Associated with these or- 
der parameters, a nonlocal unitary transformation can 
be constructed to expose explicitly the Z 2 x Z 2 symme- 
try of the Hamiltonian [y, [D, |8(. However, a nonlocal 
string order parameter that reflects correctly the hidden 
Zs+i x Zs+i topological symmetry of the higher-,!? va- 
lence bond solid has not been found [9(. 

In this paper, we introduce a novel matrix product 
state with SO(2n + 1) symmetry and show that it is 
the exact ground state of a model Hamiltonian with 
nearest neighbor interactions constructed with either the 
SO(2n + 1) projection operators or more generally the 
SU(2) spin projection operators. Unlike the valence bond 
solid state, we find that the hidden topological order in 
this class of matrix product states can be fully identi- 
fied and characterized by a set of nonlocal string order 
parameters. When n = 1, the 50(3) symmetric matrix 
product state is exactly the same as the S = 1 valence 
bond solid state and the model Hamiltonian possesses a 
hidden Z 2 X Z 2 topological symmetry [y, LD, [8(. When 



n > 1, it will be shown that the SO(2n+ 1) ground state 
possesses a hidden (Z 2 x Z 2 ) n topological symmetry. The 
breaking of this hidden symmetry leads to 4™ degenerate 
ground states with disentangled edge states in an open 
chain system. 

Let us start by considering a one dimensional lattice 
system with SO(2n + 1) symmetry. Each lattice site con- 
tains 2n + 1 basis states {\n a } ; a = 1, • • • , 2n + 1}, which 
can be rotated within the SO(2n + 1) space as follows 



L ab \n c ) = iS bc \n a )-iS ac \n b ) 1 



(1) 



where L ab (a < b) are the (2n 2 + n) generators of the 
SO(2n + 1) Lie algebra, satisfying the following commu- 
tation relations 



[L ab ,L cd ]^i(S ad L 



be 



S bc L ad 



-bd 



S bd L a 



(2) 



According to the Lie algebra, the product of any two 
SO(2n + 1) vectors can be decomposed as a sum of an 
SO(2n+l) scalar 1, an antisymmetric SO(2n + l) tensor 
2n 2 + n , and a symmetric SO(2n + 1) tensor 2n 2 + 3n : 



2n + 1 <g> 2n + 1 = 1 2n z + n ® 2rT + 3n. 



(3) 



The number above each underline is the dimension of the 
irreducible representation. 

In the spinor representation, the SO(2n+l) generators 
can be expressed as r ofc = [T a , T b ]/2i, where L Q (a = 1 ~ 
2n + 1) are the 2" x 2" matrices that satisfy the Clifford 
algebra {L a , L 6 } = 28 a b [III- For each lattice site i, if the 
following matrix state is introduced 



fh 



= E ra 



then it can be readily shown that the bond product of gi 
at any two neighboring sites have finite projection only 



in the scalar 1 and the antisymmetric 2n 2 + n subspaces 
spanned by \nf) and \nf +1 ) states, because the product 
of T a and T b can be expressed as r Q r 6 = S ab + iT ab . 
This is a special property of the SO(2n + 1) spinor rep- 
resentation constructed by Clifford algebra. By applying 
this argument to a periodic chain, we can show that the 
matrix product state defined by 



I*) 



Tr (gig 2 ■ ■ ■ 9l) ■ 



(4) 



is the exact ground state of the following SO(2n + 1) 
symmetric Hamiltonian 



H 



SO(2n+l) 



= X^2rf+3n(M+l) 



(5) 



where V 2n 2 +3n (i, j) is a projection operator that projects 
the states at sites i and j onto their SO(2n+l) symmetric 
tensor 2n 2 + in . To compute the static correlation func- 
tions of the matrix productground state ((4]), we can use 
a transfer matrix method [a, lllj . At large distance, the 
two-point correlation functions of SO(2n + 1) generators 
decay exponentially as 



(LfLf) 



exp 



|j-*l 
4 



(0) 



2w+l 
2n-3 



with the correlation length £ = 1/ In 

For the three SO(2n + 1) channels given in Eq. ([3]), 
the bond Casimir charge J2a<b(^i b + Lj b ) 2 for two ad- 
jacent sites takes the values 0, An — 2, and in + 2, 
respectively. Combining this result with the equation 
J2a<b(^i b ) 2 = ^ n an d the completeness condition of the 
projection operators, we can then express the bond pro- 
jection operator V2n 2 +3n(h j) with the SO(2n + 1) gen- 
erators as 



2n 2 +3n 



(U) 



\ * rabrab . /\ * r ab r ab\2 



a<b 



a<b 



2n+ 1 



Thus the model defined by Eq. ((5J is a. bilinear- 
biquadratic Hamiltonian in terms of the SO(2n+ 1) gen- 
erators. 

At each lattice site, the 2n+l vectors of SO(2n+l) can 
be also constructed from the S = n quantum spin states. 
In the SU(2) spin language, the last two channels in Eq. 
([3]) correspond to the total bond spin S = 1, 3, . . . , 2n— 1 
and S — 2,4, ...,2n states, respectively. Furthermore, 
it can be shown that the bond projection operators of 
SO(2n + 1) can be expressed using the spin projection 
operators Ps= m (hj) as 



T > 2n 2 +n(hj) 



m—l 



Ps=2m-l(hj), 



^2n 2 +3 



n(i,j) 



}_^ Ps=2m(i,j) 

m—l 



Thus V 2 n 2 +3n (i,j 

i and j onto the nonzero even total spin states 



is to project the spin states at sites 

Based 



on this property, we can further show that the matrix 
product wavefunction (j4|) is also the ground state of the 
following integer spin Hamiltonian 



H 






517(2) =2-^2.^ JmPs=2m(i,i+ 1) 



(7) 



with all J m > 0. This model is S'C/(2)-invariant in gen- 
eral. However, the ground state ((4]) possesses an emer- 
gent SO(2n + 1) symmetry. When all J m = 1, Hsuti) 
becomes SO(2n + l)-invariant. In this case, Hgu/ 2 ) sim- 
ply reduces to H so(2n+ i). 

It is interesting to compare Hsu (2) with the AKLT 
model of valence bond solid proposed by Affleck et. al. 



H A 



KLT 



E 



2)> 

E 



K m P s 



S, 



(8) 



i m— n+1 



with all K m > 0. The ground state of Haklt is also a 
matrix product state similar to Eq. (j4|), but gi is now 
a (S + 1) x (5 + 1) = (n + 1) x (n + 1) matrix [|]. 
These two matrix product states have different topologi- 
cal properties and belong to different topological phases 
when n > 1. Therefore Hsu(2) an d Hso(2n+i) can be 
viewed as a new family of exactly solvable quantum in- 
teger spin models to understand the internal structures 
of Haldane gap phases. 

When n = l, both Hso(2n+i) an d Hsu<2) become 
exactly the same as the S = 1 AKLT model Haklt- 
The ground state has a hidden antiferromagnetic order 
in which the up and down spins lie alternately along the 
lattice, sandwiched by arbitrary number of non-polarized 
spin states. This dilute antiferromagnetic order can be 
measured by a nonlocal string order parameter first pro- 
posed by den Nijs and Rommelsc 5], 



C"' 



lim 

-i|-M 






At: St 



s$) 



4 
9' 



(9) 



where /i = x, y or z. By performing a nonlocal unitary 
transformation [O, LD, LSI] to the spin operators with the 
following unitary operators 



U = Y[exp(inS*S*), 



(10) 



j<i 



two of the above string order parameters are converted 
into the conventional spin-spin correlation functions. The 
SU(2) symmetry of the AKLT model is then reduced to a 
discrete Z 2 x Z 2 symmetry [y, |7|, [8( . This reveals a hidden 
topological symmetry of the original model. The break- 
ing of this topological symmetry leads to the opening 
of the Haldane gap and the four-fold degenerate ground 
states in an open chain. 



Similar to the n = 1 case, the general SO(2n + 1) 
(n > 1) matrix product state (|4|) also contains interesting 
hidden antiferromagnetic orders. Since SO(2n + 1) is a 
rank-n algebra, one can always classify the states at each 
site using n quantum numbers (weights) {mi, • • • , m„} 
subjected to the constraint 



m a mp — 0, (a ^ f3). 



(11) 



Here {mi, • • • , m„} are the eigenvalues of the mutually 
commuting Cartan generators {L 12 , L 34 , . . . , t^ 2 ™- 1 . 2 "} 



-2a-l,2a| 



m a \m a ) 



(m Q = 0,±l) 



(12) 



According to Eq. (p}, ah these Cartan generators anni- 
hilate the state |?i 2 " +1 ) = |0, 0, . . . , 0) . The other basis 
states are given by 



|0- 



±1, 



■i)) = -^(\n 2 «)± l \n 2 ^)). (13) 



From the property of the Clifford algebra, the hidden 
antiferromagnetic order of the ground state \^f) can now 
be identified. In any of these m a (a = 1 ~ n) chan- 
nel, it can be shown that \m a ) is dilute antiferromag- 
netically ordered, same as for the 5 = 1 valence bond 
solid. Namely, the states of m a — 1 and m a = — 1 will 
alternate in space if all the m a = states between them 
are ignored. For example, a typical configuration of the 
ground state of the 5*0(5) system is 

mi: -0t00||000|t010t ••• 
m 2 : ■•■ t0|000ti000f0|0- 

where (t,0,i) represent \m) = (|1), |0), | — 1)) states, re- 
spectively. 

This hidden antiferromagnetic order reminds us a gen- 
eralization of the den Nijs-Rommelse nonlocal string or- 
der parameters to characterize this state. Similar to Eq. 
© of the n — 1 case [5j , the string order parameters can 
be defined as 



O ab = lim (If TT exp(i-irLf)L 

I 4 _ 4. I — krvi J- J- 



\3 



(14) 



Since the ground state is SO(2n+ 1) rotationally invari- 
ant, the above nonlocal order parameters should all be 
equal to each other. Thus to determine the value of these 
parameters, only the value of O 12 needs to be evaluated. 
In the L 12 channel, the role of the phase factor in Eq. 
(fT4"|) is to correlate the finite spin polarized states in the 
mi channel at the two ends of the string. If nonzero mi 
takes the same value at the two ends, then the phase fac- 
tor is equal to 1. On the other hand, if nonzero mi takes 
two different values at the two ends, then the phase fac- 
tor is equal to —1. Thus the value of O 12 is determined 
purely by the probability of mi = ±1 appearing at the 
two ends of the string. Since the ground state is trans- 
lation invariant, it is straightforward to show that the 



probability of the states mi = ±1 appearing at one lat- 
tice site is 2/(2n + 1) and thus O 12 = 4/(2ra + l) 2 . 

The Kennedy- Tasaki unitary transformations 1|10[) for 
n = 1 case [a, @, Q can also be generalized to ar- 
bitrary n > 1 cases. In the SO(2n + 1) Lie al- 
gebra, (L 2a - 1 < 2a ,L 2a - 1 ' 2n+1 ,L 2a ' 2n+1 ) span an 5*0(3) 
sub-algebra in which exp(i7rL 2Q,2n+1 ) plays the role of 
flipping the quantum number m Q . This exponential op- 
erator can flip the quantum numbers of m a without dis- 
turbing the quantum states in all other channels. This 
indicates that if we take the following nonlocal unitary 
transformation in the m a channel 

tt TT (■ f 2a-l,2a r 2a,2n+l\ / 1c v 

Ua = 1 1 eXP [ mL j L i )' ( 15 ) 

j<i 

then all the configurations in this channel will be ferro- 
magnetically ordered. Furthermore, by performing this 
nonlocal transformation successively in all the channels 



u=Y[u a , 



(16) 



then all the configurations of the ground state will be- 
come ferromagnetically ordered. As an example, Fig. [T] 
shows how the 50(5) matrix product state \^) is succes- 
sively changed under this nonlocal unitary transforma- 
tion. 

By applying the generalized unitary transformation 
(|16|) to the Cartan generators, it can be shown that 



ab , 



ULfU- 1 =Lfexp{iir 



i-l 



Substituting this formula to Eq. (TT4|) . we find that 

'ib _ i,- / T ab rab\ 



O ab = lim (LfLf) 



U ' 



(17) 



(18) 



Thus the nonlocal string order parameters O ab become 
the ordinary correlation functions of local operators after 
the unitary transformation. 

Under the above transformation, the symmetry of the 
original Hamiltonian -ffgo(2n+i) is reduced, and deter- 
mined by the symmetry of the unitary transformation 
operators. In the m a channel, it can be shown that 
the unitary operator U a possesses only a Z2 x Z2 sym- 
metry [y, [tJ, |8j. Therefore, the Hamiltonian after the 
transformation has a (Z2 X Z<z) n symmetry. This is the 
hidden topological symmetry of the original Hamiltonian 
Hso(2n+i), associated with the hidden topological order 
of the original matrix product state |\P). Furthermore, 
the unitary transformation (|16p breaks the translational 
symmetry. When it is applied to an open chain system, 
the hidden (Z2 x Z2Y 1 topological symmetry of the Hamil- 
tonian will be further broken, yielding 2 n free edge states 
at each end of the chain. Therefore, the open chain has 
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The exactly soluble point #mps = tan 



2n+l 



has been 
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FIG. 1: (Color online) Changes of a typical configuration of 
the 50(5) ground state under the unitary transformation de- 
fined by Eq. (|16[1 . Ui and U2 transform successively all mi 
and mi states to two diluted ferromagnetic configurations, 
respectively. 



totally 4™ degenerate ground states, which can be distin- 
guished by their edge states. 



As already mentioned, H, 



SO(2n+l) 



is a bilinear- 



biquadratic Hamiltonian in terms of the SO(2n + 1) 
generators. Actually, we can introduce a general one- 
parameter family of the SO(2n + 1) bilinear-biquadratic 
model as 



H 



E 



a<b 



E L ^+i 

\a<b 1 



(19) 

which is an extension of the quantum spin-1 bilinear- 
biquadratic model. To determine the region of the Hal- 
dane gapped phase, we need to identity several special 
integrable points. At 9\ — tan -1 -*_ , the model 111 



becomes the Uimin-Lai-Sutherland (ULS) model with an 
enhanced SU(2n + 1) symmetry, which can be solved 



by Bethe ansatz [12j. It is well-known that this model 
has gapless excitations described by SU(2n + l)i Wess- 
Zumino-Witten model [L3[. Based on the renormaliza- 
tion group approach, for < 6\, Itoi and Kato [14J found 
that the marginally relevant interaction generates the 
Haldane gap, and the transition at the ULS point be- 
longs to the universality class of the Kosterlitz-Thouless 
phase transition. 

One the other hand, using quantum inverse scattering 
methods, Reshctikhin [15] had discovered another class of 
one-dimensional quantum integrable SO(n) model, cor- 
responding to the point 9 2 — tan -1 ^n-i)' 1 > wnere there 
are also gapless excitations above the ground state. For 
n = 1, this point corresponds to the quantum spin-1 
Takhatajan-Babujian model 16], which is at the bound- 
ary between Haldane gap phase and dimerized phase. 
These rigorous results suggest that the Haldane gapped 
phase for the general model (fT9|) exists in the region 



tan 



2n-3 
(271-1)^ 



< 9 < tan~ 



1 



2ra-l 



(20) 



included. In the whole region, we expect that the system 
has an energy gap in the excitations and the ordinary cor- 
relation functions display exponentially decay. However, 
a nonvanishing string order parameter (|14[) can measure 
the breaking of the hidden topological symmetry. 

For n = 1, the spin-1 quantum antiferromagnetic 
Heisenberg model (9 = 0) is just included in this region, 
however, we find that the SO{2n + 1) Heisenberg point 
for n > 2 does not belong to the Haldane gap phase. In 
particular, when n — 2, the corresponding 50(5) antifer- 
romagnetic Heisenberg model has been used by Scalapino 
et. al. 17] to describe the 50(5) "superspin" phase on 
a ladder system of interacting electrons. Therefore, the 
ground state and low-lying excitations of the quantum 
SO(2n + 1) symmetric generalized Heisenberg model for 
n > 2 deserves further studies. 

In conclusion, we have constructed an SO(2n + 1)- 
invariant matrix product state and shown that it is the 
exact ground state of an SO(2n + l)-symmetric Hamil- 
tonian defined by Eq. ((5|) or more generally an SU(2)- 
symmetric spin Hamiltonian defined by Eq. ([7]). This 
matrix product state contains diluted antiferromagnetic 
orders in n different channels and a hidden (Zi x Z-i) n 
topological symmetry. These topological long range or- 
der can be characterized by a set of nonlocal string order 
parameters. The breaking of the [Z<x x Z<i) n topological 
symmetry leads to the opening of an excitation gap be- 
tween the ground state and the first excitation state. In 
an open chain system, the 4 n edge states become com- 
pletely disentangled and the ground states are 4™ degen- 
erate. The multiple Z 2 nature of these topological states 
suggests that they can serve as a resource of multiple 
qubits. We believe that these states, similar as for the 
5 = 1 AKLT valence bond state, can be encoded to per- 
form ideal quantum teleportation [181 ] or fault-tolerant 
quantum computation through local spin measurements. 
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